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We study the accuracy of determining the phase space quasidistribution of a single quantized 
light mode by a photon counting experiment. We derive an exact analytical formula for the error 
of the experimental outcome. This result provides an estimation for the experimental parameters, 
such as the number of events, required to determine the quasidistribution with assumed precision. 
Our analysis also shows that it is in general not possible to compensate the imperfectness of the 
photodetector in a numerical processing of the experimental data. The discussion is illustrated with 
Monte Carlo simulations of the photon counting experiment for the coherent state, the one photon 
Fock state, and the Schrodinger cat state. 
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I. INTRODUCTION 

Generation and detection of simple quantum systems 
exhibiting nonclassical features have been, over the past 
few years, a subject of numerous experimental and the- 
oretical studies. Quantum optics, provides in a natural 
way, several interesting examples of nonclassical systems 
along with practical tools to perform measurements. One 
of the simplest systems is a single quantized light mode, 
whose Wigner function has been measured in a series of 
beautiful experiments using optical homodyne tomogra- 
phy |l],H following a theoretical proposal || . 

Recently, a novel method for the determination of the 
Wigner function of a light mode has been proposed [|],|| . 
This method utilizes a direct relation between the pho- 
tocount statistics and the phase space quasidistributions 
of light. Simplicity of this relation makes it possible to 
probe independently each point of the single mode phase 
space. In a very recent experiment, an analogous method 
has been used to measure the Wigner function of the mo- 
tional state of a trapped ion || . 

In the present paper we supplement the description 
of this method with a rigorous analysis of the statisti- 
cal uncertainty, when the phase space quasidistribution 
is determined from a finite sample of measurements. The 
analysis of the reconstruction errors due, for example, to 
a finite set of recorded photons is essential in designing 
an experiment that can provide accurate reconstruction 
of the phase space quasidistribution. We will derive an- 
alytical formulae for the uncertainty of the experimen- 
tal result, which give theoretical control over the poten- 
tial sources of errors. Additional motivation for these 
studies comes from recent discussions of the possibility 
of compensating the detector losses in quantum optical 
measurements 0J^]. Our analysis will provide a detailed 
answer to the question whether such a compensation is 
possible in the newly proposed scheme. 

This paper is organized as follows. First, in Sec. || we 



briefly review properties of the quasidistribution func- 
tions that are relevant to the further parts of the pa- 
per. Then, in order to make the paper self contained, 
we present the essentials of the method in Sec. Next 
we derive and discuss the statistical error of the exper- 
imentally determined quasidistribution in Sec. |y|. The 
general results are i llus trated with several examples in 
Sec. [V[ Finally, Se 
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VI summarizes the paper. 



II. 



QUASIDISTRIBUTIONS OF A SINGLE 
LIGHT MODE 



The phase space representation of the state of a sin- 
gle quantized light mode as a function, of one complex 
variable fi, has been extensively used in quantum op- 
tics since its very beginning. Due to the noncommuta- 
tivity of the boson creation and annihilation operators, 
the Wigner function is just an example of the more gen- 
eral t-parameterized quasiprobability distributions. Such 
a one-parameter family of quasidistribution functions is 
given by the following formula M : 



W{(3;t) = \ 



d 2 £e 



m 2 /2+e p-w 



(i) 



where a and a' are the single mode photon annihila- 
tion and creation operators, and the real parameter t 
defines the corresponding phase space distribution. It 
is well known that such a parameter is associated with 
the ordering of the field bosonic operators. For example, 
t = 1,0, and —1 correspond to the normal, symmetric 
and antinormal ordering, respectively. The correspond- 
ing quasidistributions are: the P function, the Wigner 
function, and the Q function. The properties of the 
t-parameterized quasiprobability distributions are quite 
different. For example: the P function is highly singular 
for nonclassical states of light, the Wigner function is well 
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behaved for all states, but may take negative values, and 
finally the Q function is always positive definite. These 
properties reflect a general relation linking any two dif- 
ferently ordered quasidistributions via convolution with 
a Gaussian function in the complex phase space: 
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where t > t' . Thus the lower the ordering, the smoother 
the quasidistribution is, and fine details of the functions 



are not easily visible. 

This behavior can be explicitly illustrated using a su- 
perposition of two coherent states: 



(|ao) + I - «o))- 



(3) 



States of this type illustrate quantum coherence and in- 
terference between classical-like components, and are of- 
ten called quantum optical Schrodinger cats [3|JlO|]. The 
quasidistribution of this superposition is given by the for- 
mula 
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Fig. [j] shows this quasidistribution plotted for three dif- 
ferent values of the ordering parameter t. The Wigner 
function contains an oscillating component originating 
from the interference between the coherent states. This 
component is completely smeared out in the Q function, 
which can hardly be distinguished from that of a statis- 
tical mixture of two coherent states. 

Computation of a higher ordered quasidistribution 
from a given one is not a straightforward task, since the 
integral in Eq. (Q) fails to converge for t' > t. Instead, a 
relation linking the Fourier transforms of the quasidistri- 
butions can be used in analytical calculations: 
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This formula can be transformed into the form [ p~4[ 
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showing that for t < the quasidistribution is an ex- 
pectation value of a bounded operator and therefore is 
well behaved. Of course, only well behaved, nonsingular 
functions have experimental significance. 



However, its application in the processing of experimen- 
tal data would enormously amplify the statistical error 
p| . Consequently, in the case of experimentally deter- 
mined functions it is practically impossible to compute a 
higher ordered quasidistribution from the measured one. 
Therefore the ordering of the measured quasidistribution 
depends primarily on the features of a specific experimen- 
tal scheme. Optical homodyne tomography is capable of 
measuring the Wigner function, provided that the detec- 
tion is lossless. For heterodyne jl2| and double homodyne 
jL3| detection schemes the Q function is the limit, because 
in these methods an additional noise is introduced by a 
vacuum mode mixed with the detected fields. 

In our calculations a normally ordered representation 
of the quasidistribution functions will be very useful. In- 
troducing normal ordering of the creation and annihila- 
tion operators in Eq. (Q) allows to perform the integral 
explicitly, which yields: 



III. PHOTOCOUNT STATISTICS AND 
QUASIDISTRIBUTIONS 

The experimental setup proposed to determine the 
quasiprobability distribution via photon counting is pre- 
sented in Fig. ^|. The signal, characterized by an annihi- 
lation operator a is superposed by a beam splitter on a 
probe field prepared in a coherent state |a). The mea- 
surement is performed only on one of the output ports 
of the beam splitter, delivering the transmitted signal 
mixed with the reflected probe field. The quantity de- 
tected on this port is simply the photon statistics, mea- 
sured with the help of a single photodetector. We will 
consider an imperfect detector, described by a quantum 
efficiency rj. According to the quantum theory of pho- 
todetection, the probability p n of registering n counts 
by the detector is given by the expectation value of the 
following normally ordered operator: 
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where the angular brackets denote the quantum expecta- 
tion value and j7out is the operator of the time-integrated 
flux of the light incident onto the surface of the detector. 
It can be expressed in terms of the signal and probe fields 

as 
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with T being the beam splitter power transmission. The 
count statistics determined in the experiment is used to 
compute the following photon count generating function 
(PCGF): 



II(a;s) 
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where s is a real parameter, which will allow us to ma- 
nipulate to some extent the ordering of the measured 
phase space quasidistribution. We will specify the al- 
lowed range of the parameter s later. A comparison of 
the latter form of the PCGF with Eq. @ shows that it 
is directly related to a specifically ordered quasidistribu- 
tion of the signal field. After the identification of the 
parameters we obtain that the PCGF can be written as: 
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Thus the PCGF computed from the count statistics is 
proportional to the value of the signal quasiprobability 
distribution at the point j3 = — T) /Ta, determined 
by the amplitude and the phase of the probe coherent 
field |15|. We may therefore scan the whole phase space 
of the signal field by changing the parameters of the probe 
field. The main advantage of this scheme is that the mea- 
surement of the quasidistribution is performed indepen- 
dently at each point of the phase space. In contrast to 
optical homodyne tomography there is no need to gather 
first the complete set of experimental data and then pro- 
cess it numerically. 

The ordering of the measured quasidistribution de- 
pends on two quantities: the parameter s used in the 
construction of the PCGF from the count statistics, and 
the product of the beam splitter transmission T and the 
detector efficiency 77. Let us first consider the case of 
an ideal detector 77 = 1 and the transmission tending to 
one. In this limit the ordering of the quasidistribution 
approaches s. Setting s — 0, which corresponds to mul- 
tiplying the count statistics by (—1)™, yields the Wigner 
function of the signal field. However, in the limit T — > 1 
only a small fraction of the probe field is reflected to the 
detector, which results in the vanishing factor multiply- 
ing a in Eq. (O), and thus an intense probe field has to 



be used to achieve the required shift in the phase space. 
Consequently, the beam splitter should have the trans- 
mission close to one, yet reflect enough probe field to al- 
low scanning of the interesting region of the phase space. 
This decreases the ordering of the measured quasidistri- 
bution slightly below zero. A more important factor that 
lowers the ordering is the imperfectness of the detector. 
Maximum available efficiency of photodetectors is lim- 
ited by the current technology and on the level of single 
photon signals does not exceed 80% Therefore in a 
realistic setup the parameter deciding about the magni- 
tude of the product r\T is the detector efficiency rather 
then the beam splitter transmission. 

One may attempt to compensate the lowered ordering 
by setting s = 1 — 77T, which would allow to measure the 
Wigner function regardless of the value of r)T. However 
in this case the magnitude of the factor multiplying the 
count statistics in the PCGF diverges to infinity, which 
would be a source of severe problems in a real experi- 
ment. We will discuss this point in detail later, in the 
framework of a rigorous statistical analysis of the mea- 
surement. 

Let us close this section by presenting several exam- 
ples of the photocount statistics for different quantum 
states of the signal probed by a coherent source of light. 
The most straightforward case is when a coherent state 
|ao) enters through the signal port of the beam splitter. 
Then the statistics of the registered counts is given by 
the Poisson distribution: 
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where J{cx.q) = r]T\/3— ao\ 2 is the average number of reg- 
istered photons. When the measurement is performed at 
the point where the quasidistribution of the signal field 
is centered, i.e., (3 — ao, the fields interfere destructively 
and no photons are detected. In general, for an arbi- 
trary phase space point, the average number of registered 
photons is proportional to the squared distance from «o- 
Averaging Eq. (12) over an appropriate P representation 
yields the photocount statistics for a thermal signal state 
characterized by an average photon number n: 
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where L n denotes the nth Laguerre polynomial. 

A more interesting case is when the signal field is in 
a nonclassical state. Then the interference between the 
signal and the probe fields cannot be described within 
the classical theory of radiation. We will consider two 
nonclassical states: the one photon Fock state and the 
Schrodinger cat state. The count statistics can be ob- 
tained by calculating the quantum expectation value of 
PCGF over the considered state and then expanding it 
into the powers of (s + l)/(s — 1). 
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The photocount distribution for the one photon Fock 
state |1) can be written as an average of two terms with 
the weights r]T and 1 — rjT: 

p W = ^(n-J(0)) 2 »^ e -^) 

+(1 _ rfD^MLe-m . (i4) 

The second term corresponds to the detection of the vac- 
uum signal field. Its presence is a result of the detector 
imperfectness and the leakage of the signal field through 



the unused output port of the beam splitter. This term 
vanishes in the limit of rfT — > 1, where the Wigner func- 
tion is measured in the setup. The first term describes 
the detection of the one photon Fock state. In Fig. |^(a) 
we show the statistics generated by this term for different 
values of (3. If the amplitude of the probe field is zero, 
we detect the undisturbed signal field and the statistics 
is nonzero only for n = 1. The distribution becomes flat- 
ter with increasing (3. Its characteristic feature is that it 
vanishes around n w J(0). 

For the Schrodinger cat state defined in Eq. (|J) the 
photocount statistics is a sum of three terms: 
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The first two terms describe the two coherent compo- 
nents of the cat state, whereas the last one contributes 
to the quantum interference structure. In Fig. ||(b) we 
plot the photocount statistics for different values of (3 
probing this structure, in the limit r\T — > 1. The four 
values of (3 correspond to the cosine function in Eq. (Q) 
equal to 1, 0, —1, and 0, respectively, for t = 0. It is seen 
that the form of the statistics changes very rapidly with 
(3. This behavior becomes clear if we rewrite the PCGF 
(Eq. (|l0|)) for s = in the form 
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1=0 



1=0 



showing that in order to obtain a large positive (negative) 
value of the quasidistribution the photocount statistics 
has to be concentrated in even (odd) values of n. 



IV. STATISTICAL ERROR 

The proposed measurement scheme is based on the re- 
lation between the quasidistributions and the photocount 
statistics. In a real experiment the statistics of the de- 
tector counts cannot be known with perfect accuracy, as 
it is determined from a finite sample of N measurements. 
This statistical uncertainty affects the experimental value 
of the quasidistribution. Theoretical analysis of the sta- 
tistical error is important for two reasons. First, we need 
an estimation for the number of the measurements re- 
quired to determine the quasidistribution with a given 
accuracy. Secondly, we have seen that one may attempt 
to compensate the imperfectness of the detector and the 
nonunit transmission of the beam splitter by manipulat- 
ing the parameter s. However, in this case the magnitude 



of the factor multiplying the count statistics diverges to 
infinity. This amplifies the contribution of the tail of 
the distribution, which is usually determined from a very 
small sample of data, and may therefore result in a huge 
statistical error of the final result. Our calculations will 
provide a quantitative analysis of this problem. 

Due to extreme simplicity of Eq. ( [io| ) , linking the count 
statistics with the quasidistributions, it is possible to per- 
form a rigorous analysis of the statistical error and obtain 
an exact expression for the uncertainty of the final result. 
We will assume that the maximum number of the pho- 
tons that can be detected in a single measurement cannot 
exceed a certain cut-off parameter K. Let us denote by 
k n the number of measurements when n photons have 
been detected, n = 0, 1, . . . , K. The set of fc n 's obeys the 
multinomial distribution 



V{koM 



,k K ) = 



Nl 



x Po ^ 1 



k lk 1 l...[N-j:^ k n )l 

/ k x^-£to fc " 

V n=0 ) 



(17) 



The measured count statistics is converted into the ex- 
perimental PCGF: 
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which is an approximation of the series defined in 
Eq. (|l(]). In order to see how well n exp approximates 
the ideal quantity we will find its mean value and its 
variance averaged with respect to the distribution (0). 
This task is quite easy, since the only expressions we need 
in the calculations are the following moments: 
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We use the bar to denote the statistical average with re- 
spect to the distribution V(ko, . . . kx). Given this result, 
it is straightforward to obtain: 
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and 
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The error introduced by the cut-off of the photocount 
statistics can be estimated by 
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The variance <5n 2 xp is a difference of two terms. The sec- 
ond one is simply the squared average of n cxp . The first 
term is a sum over the count statistics multiplied by the 
powers of a positive factor ((s + l)(s— l)) 2 . If s > 0, this 
factor is greater than one and the sum may be arbitrarily 
large. In the case when the contribution from the cut tail 
of the statistics is negligible, i.e., if AT ^ oo, it can be 
estimated by the average number of registered photons: 
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Thus, the variance grows unlimited as we probe phase 
space points far from the area where the quasidistribu- 
tion is localized. Several examples in the next section 
will demonstrate that the variance usually explodes much 
more rapidly, exponentially rather than linearly. This 
makes the compensation of the detector inefficiency a 
very subtle matter. It can be successful only for very 
restricted regions of the phase space, where the count 
statistics is concentrated for a small number of counts 
and vanishes sufficiently quickly for larger n's. 

Therefore, in order to ensure that the statistical error 
remains bounded over the whole phase space, we have to 
impose the condition s < 0. Since we are interested in 
achieving the highest possible ordering of the measured 
quasidistribution, we should consequently set s ~ 0. For 
this particular value the estimations for the uncertainty 
of n oxp take a much simpler form. The error caused by 



the cut-off of the count distribution can be estimated by 
the "lacking" part of the probability: 
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which shows that the cut-off is unimportant as long as 
the probability of registering more than K photons is 
negligible. The variance of n exp is given by 
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Thus, the statistical uncertainty of the measured qua- 
sidistribution can be simply estimated as 1/yN multi- 
plied by the proportionality constant given in Eq. (|ll"|). It 
is also seen that the uncertainty is smaller for the phase 
space points where the magnitude of the quasidistribu- 
tion is large. 

Analysis of the statistical error for the recently demon- 
strated measurement of the Wigner function of a trapped 
ion is different, as the analog of p n is not determined by a 
counting-like method (l^] . A quantity that is detected in 
the experimental setup is the time-dependence of the flu- 
orescence light produced by driving a selected transition. 
This signal is a linear combination of the populations of 
the trap levels with known time-dependent coefficients, 
and the p n 's are extracted from these data by solving the 
resulting overdetermined set of linear equations. Conse- 
quently, their statistical error is not described by a simple 
analytical formula. Additional source of error is the un- 
certainty of the phase space displacement, which is per- 
formed by applying an rf field. An element of our analysis 
that can be directly transferred to the case of a trapped 
ion is the effect of the cut-off of the count statistics. 



V. EXAMPLES 

We will now consider several examples of the recon- 
struction of the quasidistributions from the data collected 
in a photon counting experiment. Our discussion will be 
based on Monte Carlo simulations compared with the 
analytical results obtained in the previous section. 

First, let us note that the huge statistical error is not 
the only problem in compensating the detector ineffi- 
ciency. If s > 0, the sum (|2^) does not even have to con- 
verge in the limit of K — > oo. An example of this patho- 
logical behavior is provided by a thermal state, which has 
been calculated in Eq. (|T3|). For the zero probe field we 
obtain 
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which shows that if s > 0, then for a sufficiently in- 
tense thermal state the magnitude of the summand is 
larger than one and consequently the sum diverges, when 
K — > oo. This behavior is due to the very slowly van- 
ishing count distribution and it does not appear for the 



other examples of the count statistics derived in Sec. |TJ. 

In Fig. U we plot the reconstructed quasidistributions 
for the coherent state \ctQ = 1) and the one photon Fock 
state. Due to the symmetry of these states, it is sufficient 
to discuss the behavior of the reconstructed quasidistri- 
bution on the real axis of the phase space. The cut-off 
parameter is set high enough to make the contribution 
from the cut tail of the statistics negligibly small. The 
quasidistributions are determined at each phase space 
point from the Monte Carlo simulations of N = 1000 
events. The grey areas denote the statistical uncertainty 
calculated according to Eq. (^l|). The two top graphs 
show the reconstruction of the Wigner function in the 
ideal case rjT =1. It is seen that the statistical error is 
smaller, where the magnitude of the Wigner function is 
large. In the outer regions it approaches its maximum 
value 1/VjV. The effect of the nonunit rjT is shown in 
the center graphs. The measured quasidistributions be- 
come wider and the negative dip in the case of the Fock 
state is shallower. In the bottom graphs we depict the 
result of compensating the nonunit value of rjT by set- 
ting s = 1 — rjT. The compensation works quite well in 
the central region, where the average number of detected 
photons is small, but outside of this region the statisti- 
cal error explodes exponentially. Of course, the statisti- 
cal error can be suppressed by increasing the number of 
measurements. However, this is not a practical method, 
since the statistical error decreases with the size of the 
sample only as 1/ \fN. 

The reconstruction of the interference structure of the 
Schrodinger cat state is plotted in Fig. [|. This struc- 
ture is very fragile, and its precise measurement requires 
a large sample of events. In the case of the presented 
plot, N — 5000 simulations were performed at each phase 
space point. Comparison of the top and the center graphs 
shows how even a small imperfectness in the detection de- 
stroys the interference pattern. The data collected in an 
nonideal setup can be processed to recover the Wigner 
function, but at the cost of a significantly larger statisti- 
cal error, as it is shown in the bottom graph. Outside the 
interference structure, we again observe the exponential 
explosion of the dispersion due to the increasing intensity 
of the detected light. 

Finally, let us look at the effect of cutting the statistics 
at a finite value. Fig. || shows the Wigner functions for 
the one photon coherent and Fock states reconstructed 
from the count distributions cut at K = 11. We per- 
formed a large number of N = 10 4 simulations in order 
to get the pure effect of the cut-off that is not spoiled 
by the statistical uncertainty. The grey areas show the 
cut-off error, estimated using Eq. (p4|). The reconstruc- 
tion works well as long as the probability of detecting 
more than K photons is negligible. When the average 



number of incident photons starts to be comparable with 
the cut-off parameter, "ghost" peaks appear. When we 
go even further, the Wigner function again tends to zero, 
but this is a purely artificial effect due to the virtually 
vanishing count distribution below K. 



VI. CONCLUSIONS 

The newly proposed scheme for measuring the 
quasiprobability distributions is based on a much more 
direct link with the data collected in an experiment than 
optical homodyne tomography. Furthermore, each point 
of the phase space can be probed independently. The 
simplicity of the underlying idea allowed us to discuss 
rigorously the experimental issues, such us the statisti- 
cal uncertainty and the effect of the cut-off of the count 
statistics. 

We have seen that it is, in general, not possible to 
compensate the detector inefficiency. This conclusion 
may seem to be surprising, since it was shown that the 
true photon distribution can be reconstructed from pho- 
tocount statistics of an imperfect detector, provided that 
its efficiency is larger that 50% [||. Moreover, a simple 
calculation shows that an application of this recipe for 
recovering the photon distribution is equivalent to set- 
ting the parameter s above zero. The solution of this 
seeming contradiction is easy: the proposed recipe can 
be safely used as long as we are interested in the pho- 
ton statistics itself. However, the statistical errors of the 
reconstructed probabilities are not independent, and an 
attempt to utilize these probabilities to calculate another 
physical quantity may lead to accumulation of the sta- 
tistical fluctuations and generate a useless result. The 
simple alternating series ^2 n {~l) n p n is an example of 
such a quantity. 
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FIG. 1. Quasidistributions representing the Schrodinger cat state for cto 
t = 0, -0.1, and -1. 
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FIG. 2. The experimental setup proposed to measure quasidistribution functions. The signal a is superposed by the beam 
splitter BS on a coherent state \a). The photon statistics of this superposition is measured by the photodetector PD. 
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(a) one photon Fock state (b) Schrodinger cat state 
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FIG. 3. The photocount statistics of (a) the one photon Fock state and (b) the Schrodinger cat state for e*o = 3i, shown for 
several values of the rescaled probe field amplitude (3 — — T)/Ta in the limit rfT — 1. 



8 



Coherent state 



Fock state 






/ r/T = 1 
' s = 




77T = 0.8 
s = 






-i 


/ 

' s = 0.2 



-4-3-2-10 1 2 3 4 

/? 

FIG. 4. Reconstruction of the quasiprobability distributions of the coherent state |qo = 1) (left) and the one photon Fock 
state (right) from N — 1000 events. The solid lines are the analytical quasidistributions and the grey areas mark the statistical 
dispersion. 
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FIG. 5. Reconstruction of the interference structure of the Schrodinger cat state for ao = Si from N = 5000 events at each 
point. 
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